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ABSTRACT 

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion 
pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with 
photometric redshifts, spanning between z = 0.45 and z = 0.65, constructed from the SDSS using 
methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a 
volume of 3/i~^Gpc^, making it the largest volume ever used for galaxy clustering measurements. We 
describe in detail the construction of the survey window function and various systematics affecting our 
measurement. With such a large volume, high precision cosmological constraints can be obtained given 
a careful control and understanding of the observational systematics. We present a novel treatment of 
the observational systematics and its applications to the clustering signals from the data set. In this 
paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with 
an accuracy of ~ 15% with bin size of 6i = 10 on scales of the Baryon Acoustic Oscillations (BAO) 
(at £ 40 — 400 ). We also apply corrections to the power-spectra due to systematics, and derive 
cosmological constraints using the full-shape of the power-spectra. For a flat ACDM model, when 
combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) 
and Hq constraints from using 600 Cepheids observed by Wide Feild Camera 3 (WFC3) (HST) , we find 
Qa = 0.73 ± 0.019 and Hq to be 70.5 ± 1.6 s-^Mpc"^km. For an open ACDM model, when combined 
with WMAP7 + HST, we find = 0.0035 ± 0.0054, improved over WMAP7+HST alone by 40%. 
For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 ± 0.078, and Hq 
to be 71.3 ± 1.7 s~^Mpc~^km, which is competitive with the latest large scale structure constraints 
from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 
2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). We also find that systematic-corrected 
power-spectra gives consistent constraints on cosmological models when compared with pre-systematic 
correction power-spectra in the angular scales of interest. The SDSS-III Data Release 8 (SDSS-III 
DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, 
cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. Here, we refer 
to our companion papers (Seo et al. 2011, de Putter et al. 2011) for further investigations using the 
clustering data. Our calculation of survey selection function, systematics maps, likelihood function 
for COSMOMC package will be released at |http: / /portal. nersc.gov / project /boss / galaxy / photoz /] 
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1. INTRODUCTION 

The distribution of light in the Universe has long been 
used as a probe into the structure of the Universe. Ein- 
stein wrote of the distribution of stars as possibly being 
uniform on average over large enough distances in 1917 
when he discussed the structure of the Universe. Hubble 
tested the uniformity of distribution of faint nebulae in 
1926. As the structure of the Universe unfolds, distribu- 
tion of light from objects such as galaxies has remained a 
powerful cos mological probe (Peebles 1973; Groth 1973j 
[Wang et al."T999; Hu 19991 [Eisenstein et al. 1999|). 

Smoothed over large scales, we expect galaxy density 
to have a simple relationship to the underlying matter 
density; this implies that the clustering of galaxies at 
large scales is directly related to the clustering of the un- 
derlying matter and is thus a sensitive probe of both the 
initial conditions of the Universe and its subsequent evo- 
lution. It is therefore not surprising that a large fraction 
of the effort in observational cosmology had been devoted 
to measuring the spatial distribution of ga laxies, as in 

,The APM 

^ ^ ^The DEEP 

^ (DEEP, |Koo||1998D VIMOS-VLT'Deep Survey 
J), Two-Decree Fie ld Galaxy 



to measuring tne spatiai QisinpuTiion oi ^a ia: 
the CfA Redshift Surv ey (iHuchra et al. | [T983 | ) , 
GalaxySurvey (APM JiVladdox et al.| |l990|, T 



(VVD57|'Le Fevre"etar 



Redshift Survey (2dl^'GRS, [Cole et aL 2005), The Two 
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Micron Ah sky Survey (2MASS, Skrutskie et a L||2006D , 
COSMOS dScoville et aL 2007), Can ada- France- Hawai i 
Telescope Legacy Surve3t^(CFHTLS, |llbert et al.j2006|), 
Galaxy And Mass Assembly survey (GAMA, 'Driver et 



al. 2009), The WiggleZ Survey (Blake et al. 20WlWe 
Sloan Digital Sky Survey (SDSS, [York et al.||20TO| . By 



2008, the SDS^has probe d ^1.5 GpC" with gala xies 



while the current SDSS-III( |Eisenstein et al.||201lD wih 
finish surveying ~ 15 Gpc^ in 2014. The planned Large 
Synoptic Survey Telescope (LSSTj^wih observe - 1000 
Gpc^ of the Universe. 

Hidden in the ever-increasing volume of surveyed Uni- 
verse, is the wealth of cosmological information that had 
not been fully exploited. In particular, the large scale 
clustering of any mass tracer, usually characterized by 
its power-spectrum, in the Universe contains three fea- 
tures that are of significant interest to contemporary cos- 
mologists. The first distinguishing feature is oscillations 
in the power- spectrum caused by acoustic waves in the 
baryon-photon plasma before hydrogen recombination at 
1000, called Ba r yon Acoustic Oscillations (hereafter 



BAO) (.Peebles fc Yu^ 1MQ| [Sunyaev fc Zeld^^^^ 



Bond <^ ljt'statiiiou||l0841 |i^Qltzman||1080[ |Hu fc Wlntel 
1996[ [Eisenstein Hu||1998| ). The BAO technique has 
emerged as the new precision cosmology probe, espe- 
cially in discerning the properties of this unknown dark 
component of the Universe "Dark Energy". The BAO 
was first observed in early 2005 both in the SDSS Lu- 
minous Red Ga laxy sample ([Eisenstein et al.||2005 ), the 
2dFGRS data (|Cole et al.|[2O05l) and in 2006 by using 



photometric Luminous Galaxies (LRGs) in 3500 deg of 
SDSS (Pad manabhan et al. 2007) . However, neither of 
the signals were strong enough to place strong cosmologi- 
cal constraints via BAO. Second, the largest scales of the 
power-spectrum can be used to constrain the primordial 
potenti al of the Un iverse, thus testing infiation. In par- 
ticular, iDalaTeFal. (2008) has pointed out the relation- 
ship between non-gaussianity of the potential in the early 
Universe (due to various possible infiationary scenarios) 
and the large scale power of mass tracer in the Universe. 
Finally, at /c ~ 0.01/iMpc~^, the power spectrum turns 
over from a /c^ slope (for a scale invariant spectrum of 
initial fiuctuations), to a spectrum, caused by modes 
that entered the horizon during radiation-dominated era 
and were therefore suppressed. The precise position of 
this turnover is thus determined by the size of the horizon 
at matter-radiation equality. It corresponds to a physi- 
cal scale determined by the total matter (I^m^^) densi- 
ties and radiation densities {Q^h'^). In particular, with a 
large survey such as SDSS, various groups had used the 
large scale power- spectrum to put stringent constraints 
on cosmologi cal par ameters, most notably Zehavi et al j 
(2002); Te^ma rk et al.| (|2Q04|); Eisenstein et al. (2005); 
[Padmanabhan et al. "(2007); rercival et al., ( |2010| ); (Reid| 



http:/ /www . cfht.hawaii.edu/Science/ 
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http://www.lsst.org 
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et al. (2010). 



The SDSS has now surveyed 14, 555 deg , and with ap- 
propriate photometric selection, we can construct a large 
unifor m sample of the photometric luminous red galaxies 
(Ross et al. 2011), and their photometric redshifts that 
can be easily calibrated using the acquired spectroscopic 
redshifts of a uniform sub-sample 10%) of the pho- 
tometric galaxies. This approach allows the possibility 
of using both standard rulers (from the turn over scale 
of power spectrum, and also the baryon acoustic oscilla- 
tions) to acquire cosmological constraints. 

We make use of this opportunity to derive one of 
the most accurate measurements of the galaxy angu- 
lar power-spectra achieved to date. We start with the 
five b a nd ima^in^ of the SDSS- III DR8; |Aihara et al.| 



(2011); Eisenstein et al. (|2011l)), and photometrically 



select a sample of luminous red galaxi es, follow ing the 
CMASS galaxy selection detailed in White et al.| ( |2Qllp ; 
the details of the construction of the sample and the red- 



shift distribution is described in Ross et al. (2011). We 



then measure the angular clustermg power spectra as a 
function of redshift with an optimal quadratic estimator, 
which is proved to provide the best statistical error-bar 
when the field is Gaussian. The galaxy density field is 
not Gaussian on small scales, due to non-linear evolution; 
however, at relatively large scales, which are the scales 
we are concerned here, the field is close to Gaussian. We 
will discuss this issue in detail in the paper. With such a 
large volume of data, we realize that the effects of large 
scale systematics are not negligible. To gauge and cor- 
rect the effects of large scale systematics, we develop a 
novel method in correcting the large scale systematics 
given that we know the list of possible systematics. We 
construct the maps of various systematics, and calcu- 
late their cross-correlation with the galaxy density, the 
systematic auto-correlations and cross-correlations. We 
can then correct for these systematics applying this new 
method. 

The paper is organized as follows : Sec. |2] describes 
the construction of the sample; Sec. [3] then presents the 
theory and measurement of the angular power spectra; 
Sec. [4] discusses the various potential systematics and 
the novel method applie d in correcting for the observa- 
tional systematics. Sec. |6.1| describes the validation of 



the cosmological parameter fitting method, and Sec. 
summarizes the cosmological constraints themselves. We 
conclude in Sec. [H 

Wherev er not explicitly me ntioned, we assume best fit 
WMAP7 (Larson et al. 201ip parameters (from model 
Icdm+sz+lens as specified on LAMBDA website rj), 
which are: a fiat ACDM cosmology with Qm = O.zUo, 
= 0.0449, h = 0.71, and ag = 0.801. 

2. THE DATA 
2.1. Observations 



(2000); Eisenstein et al. (2011)) mapped over a quarter 



http:/ /lambda. gsfc.nasa.gov/product/map/dr4/parameters.cfm 



of the sky using the dedicated Sloan Foundation 2.5 m 
telescope located at Apache Point Observatory in New 



Mexico (Gunn et al. |2QQ6|). A drift-scanning mosaic 



CCD camera ( |Gunn et al. [ |l998l [Gumi et al.||2006|) im- 



aged th^_sk\M^^ve[^Hoto band-passes ( [l^ukugit a 



et al. 1996 Smith et al. 2002) to a limiting magnitude 
of rc^ 22.5. 'I'he imaging data were processed through 
a series of pipelines that perform astromet ric calibration 

Pier et al. ^^ 2003^ , photometric red uction ([Lupton et al. 

100 lb and photometric calibration (iFadmanabhan et al. 



2008). In particular, Baryon Oscillations Spectroscopic 



SHRTey (hereaf ter BOSS) which is a part o f SDSS III 
(Eisenstein et al.||2011[ [Aihara et al. 2011), has com- 
pleted an additional 3000 square degrees of imaging and 
is now obtaining spectra of a selected subset of 1.5 mil- 
lion galaxies. The targets are assigned to spectroscopic 
plates (tile s) u s in^ a n adaptive tiling algorithm based on 
[Blanton et aT] ( |2003p , and observed with a pair of fiber- 
lecl spectrograpns! 

The availability of large uniform photometric data-set 
prompted the start of this project, thus a series of pa- 
pers, starting wi th the generation of the photometric red- 
shifts catalog by 'Ross et al. (2011), which uses 112,778 of 



BOSS spectra as a trainmg sample for photometric cat- 
alog. The photometric redshift catalog contains over 1.6 
million objects, and 900,000 of these objects lie within 
our imaging mask and the selected redshift range (0.45 
< z < 0.65). The redshift range is selected so that it 
is nearly completely independent from DR7 analysis of 
LRG c lustering using spectr oscopy which stops at z 0.4 
dReid et al. 2010 ; Percival et al.| |2010D; this allows for 
possibility trivial comBination oT likelihoods. These 
galaxies are among the most luminous galaxies in the 
universe and trace a large cosmological volume while hav- 
ing high enough number density to ensure shot-noise is 
not a dominant contributor to the clustering variance. 
The majority of the galaxies have spectral energy distri- 
butions (~ 85% , see Masters et al.^ (2011) and private 
communication with the B(J!S!S galaxy-evolution group) 
that are distinctive of old stellar populations. 

2.2. Defining Luminous Red Galaxies 
We make use of the CMASS sa mple from BOSS , which 



is defined in [ White et al^ ( |201lD and [Ross et al.| ( |201lD 
and we write down the criteria here again tor conve- 
nience: 

17.5 < icmod < 19.9 

^mod '^mod ^ 2 

> 0.55 

ifiber2 < 21.7 

^cmod < 19-86 - 
cii > 1.6 



-1.6 X {d^ -0.8) 



The Sloan Digital Sky Sur vey (SDSS; [York et al] where 



d_\_ = {r — i) — {g — r)/8 ^ r — i 

C|| = 0.7 * (^ - r) + 1.2 * (r - i - 0.1^ 



(1) 



(2) 
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The magnitude s denoted by "cm od" are "cmodel mag- 



nitudes" (see White et al. (2011) for more discussions), 
and the colors are defined with model magnitudes, ex- 
cept for Zfiber2, which is the m agnitud e in the 2 spec- 
troscopic fiber (| St ought on et al. 2002; |Abazajian et al. 
|2004). Note that we applied ifiber2 < 21.7, although the 
current BOSS target selection has moved the limit from 
21.7 to 21.5. All magnitud e s are extinction corrected 
using maps of Schlegel et al. (1998). 

In addition to constructmg galaxy density maps, we 
created several additional maps that we use to reject 
regions heavily affected by sample systematics such as 
poor sky or stellar density, and to make sure our final 
power-spectra are free of systematics. These include (i) 
a map of the full width at half-maximum (FWHM) of 
the point-spread function (PSF) in r band; (ii) a map 
of stellar density (18.0 < r < 18.5 stars); (iii) a map of 
sky brightness in z-band in nanomaggiea^j/arcsec^ ; (iv) 
3 map of the color off sets in u — g — r and r — i from 
Schlafiy et al] ( |2QTq1 ); (v) a map of Galactic extinction 



simply resca led from, the extinction maps from |Schlegel| 



et al. (1998). 



2.3. Angular and Redshift Distributions 

To interpret the clustering of any sample, one must 
characterize the expected distribution of the sample as 
if it is completely random. This involves understanding 
both the angular and radial selection function in addition 
to the expected galaxy density, which is characterized by 
its mean density. 

To characterize the angular window function, we gen- 
erate the complete angular mask of the survey following 
the procedures described below. The observed sky is de- 
fined as a union of all fields. Determining the window 
function requires identifying the fields that cover each 
position on the sky and deciding which of those fields 
should be considered primary at that position. There is 
a unique set of disjoint polygons on the sky defined by all 
the field boundaries, which are calculated us ing MAN- 
GLE package F} Hamilton] |l993j |Hamilton fc Tegmark 
2004 Swanson et al.,,2008) and each held can be divided 
mto multiple polygons. We now must decide which fields 
are primary for each polygon in the sky; the process is 
described in|Aihara et al. (20TT| in detail. Once we de- 
termine which fields are primary for all the polygon in 
the sky, we make a cut on the field observing conditions 
(SCORE >= 0.6; for more details on SCORE, see |Aihara 
et al.|(|201l|) or the SDSS-III webpageFI. We now have 



an unified MANGLE polygon file that mcludes all the 
fields that are imaged in the entire SDSS footprint, with 
the correctly assigned primary fields with good observ- 
ing conditions. We call this the "full imaging mask", as 
plotted in Figure [T] The color in Figure [T] represents the 
date at which the imaging data was taken. The striped 
pattern perpendicular to the scanning direction is eas- 

http:/ /data.sdss3.org/datamodel/glossary.html#nanomaggies 

http:/ /spae. mit.edu /molly /mangle 

sdss3.org 




Fig. 1 . — The full imaging angular mask in equatorial coordinate 
system after generating a unique set of all polygons that contains 
primary fields with good observing conditions. The colors repre- 
sents the Modified Julian Date of observation of each field. 

Angular Mask 




Fig. 2. — The preliminary imaging mask after applying primary 
selection cuts such as cuts on seeing and the bright star mask on 
the full imaging angular mask. 

ily visible, and we can also see that the north and south 
Galactic caps are observed at significantly different epoch 
of the survey. This provides hint as to what potential ob- 
servational systematic effects would look like. To create 
a more restrictive mask which is catered towards photo- 
metric red galaxies, we proceed to exclude re gions where 
E(B-V) > 0.08 (Scranton et al. 2002; Ross et al. | |2QQ6l 
[Padmanabhan et al.||2007[ |Ho et al.||2008| ), which Ulil- 
most identital to > 0.2, when seeing m the z-band ex- 
ceeds 2.0^' in FWHM, and maskin g refflons aroun d stars 
in the Tycho astrometric catalog ( H0g et aL]|200 0). The 
final angular selection function covers a solid angle of 



11,000 square degrees, and is shown in F ig. [2| 
Applying the selection criteria in Secl2.2|to th€ 



e 14,555 

deg^ of photometric SDSS imaging considered in this pa- 
per yields a catalog of approximately 1,500,000 galax- 
ies. Applying the angular selection function as shown 
in Fig. [2] to the Ross et al. (2011) photometric redshift 
catalog yields a sample of 872,921 objects, 96% of which 
are believed to be galaxies (3% are stars, and 1% are 
quasars, according to statistics gather ed in the spectro- 
scopic sub-sample; Ross et al. (2011)). For every ob- 
ject, photometric redshift s and probabilities of being a 
galaxies were determined using the ANNz Neura l Net- 
work (Col lister fc Lahav| |2004; Firt h et al.p 003). The 
calibrati on and accu r acy y these data are discussed in 
detail in Ross et al. ( 2011[ ). In the range considered in 
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Fig. 3. — The redshift distribution of the photometric CMASS 
sample when we match the objects with an unbiased sub sample 
from SDSS-III BOSS. 




Fig. 4. — The photometric vs spectroscopic redshift distribution 
of 112,778 of SDSS-III BOSS CMASS galaxies. 

this paper, the redshifts have cahbrated errors ~ 0.04 
at z ~ 0.45 and ~ 0.06 at z ~ 0.65. We pixehze 
these galaxies as a weighted (with the probabihties of 
being a galaxy) numbe r over density. Sg = Sn/n^ onto a 
HEALPix pixelization ( [Gorski et al |l99"9| ) of the sphere, 
with 12,582,912 pixels over the whole sphere (HEALPix 
resolution 10, nside=1024), each pixel covers a solid an- 
gle of 11.8 arcmin^. These pixelized maps are used di- 
rectly to compute the angular power-spectra using opti- 
mal quadratic estimator. The optimal quadratic estima- 
tor does not down-sample input pixelized maps, rather, 
it computes the covariance matrix directly from these 
pixelized maps, and this will be discussed further in |3.4| 
The sample is divided into 6 photometric redshift slices 
of thickness Az = 0.05 starting at z = 0.4 for CMASS 
sample (CMASS through CMASS 5, see Table[l]for de- 
tails), and the underlying redshift distributions lor each 
slice are calculated using BOSS spectroscopic redshift of 
the same sample. The redshift distribution of the sam- 
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0.525 
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0.575 
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CMASS 4 
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TABLE 1 

Descriptions of the 6 = 0.05 redshift slices; Zmid is 

THE MIDPOINT OF THE REDSHIFT INTERVAL. BlAS PARAMETERS 

ARE DEDUCED FROM MARGINALIZING OVER ALL THE OTHER 
COSMOLOGICAL PARAMETERS (aND A FREE SHOT NOISE TERm) 
FROM COMBINING WMAP 7 + HST + DR8 ANGULAR 
POWER-SPECTRA LIKELIHOOD USING ONLY 30 < ^ < 150 

MULTiPOLES. The first and last bins are dropped from 

HERE ON DUE TO THE SMALL NUMBER OF GALAXIES IN THOSE 
BINS. 

pie is plotted in Fig. [3] We can see that although the 
majority of the objects in one photometric-redshift bin 
is in their corresponding true redshift bin, a significant 
fraction of them fall into neighboring bins. The compar- 
isons of these photometric redshifts to the spectroscopic 
redshifts (obtained via SDSS III spectra) are plotted in 
Fig. |4j while properties of the different slices are summa- 
rized in Table [T] We see that the numbers of galaxies in 
both the first and the last bins are significantly smaller 
than the others, therefore, we decide to drop these 2 bins 
from our analysis. This decision is also facilitated by the 
fact that we wish to have a nearly independent s ample 
from the i Reid et al.l ( |2010[ ); [Percival et al] ( |2010| ) LRG 
clustering analysis, thus allowing tor simple combination 
of their likelihoods in the cosmological parameter analy- 
sis. 

2.4. Sample Systematics 

There are a number of potential systematic effects in 
photometric samples that contaminate clustering: stel- 
lar contamination and obsuration, seeing variations, sky 
brightness variations, extinction and color o ffsets (such 
as tho se described in'Schlafiy et al. ( |2010| )). [Ross et aL 
(|2011|) had extensive discussions on these potential sys- 
tematics; we will concentrate on the particular effects 
from various systematics on the angular power-spectra 
in the range of scales that affects our science analysis. 

The above cuts remove only parts of the sky that are 
significantly affected by extinction and seeing variations. 
With such a large sky coverage, an accurate determina- 
tion of the angular power-spectra of the the large scale 
tracer is only possible through an thorough understand- 
ing of the systematics. However, if we only retain parts of 
the sky that have the minimum systematics effects; we 
must remove most of t he co verage, as we have demon- 
strated in Ross et al.| (|2011|). Therefore, we developed 
a novel way of dealing with residual sample systematics 
which we will discuss in Sec. [5] in detail. 

3. THE ANGULAR POWER SPECTRUM 

As was noted in the introduction, the angular power- 
spectrum contains information of both the growth and 
the expansion of the Universe through two standard 
rulers of the Universe: the baryon acoustic oscillations 
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and the matter-radiation equality turn-over scale; thus 
the shape of the power- spectrum. In this section, we will 
describe both the theory and the computation of angular 
power spectrum. 

3.1. From galaxy distributions to angular 
power- spectrum 

The intrinsic angular galaxy fluctuations are given by: 



9(0) 



dzb{z)N{z)S{x{z)e,z), 



(3) 



where b{z) is an assumed scale-independent bias factor 
relating the galaxy overdensity to the mass overdensity, 
i.e. Sg = bS^ ^{^) is the normalized selection function, 
and x(^) is the comoving distance to redshift z. We focus 
on the auto power-spectrum of the galaxies: 



C^'W = 11 k^dkP{k)[g]e{k)[g],_ 



ik) 



(4) 



where P{k) = P(/c, z = 0) is the matter power spectrum 
today as a function of the wave number /c, and the func- 
tion [g]i is 

[g], {k) = J dzb,{z)N{z)D{zMkx{z)) (5) 

The Limber approximation, which is quite accu- 
rate when i is not too small {£ > 10), can be 
obtained from Eq. Q by setting P{k) = P{k = 
{i + l/2)/x{z)) and usmg the asymptotic formula that 

(2/^) / edkn{kx)Mkx') = {1/X')S{X - X') (when I » 
1). We find that the substitution = + 1/2) /x{z) 
is a better approximation to the exact expressions than 
k = i/x{z)' Note that ji{x) is the order spherical 
Bessel function. On large scales where the mass fiuctu- 
ation (5 <C 1, the perturbations grow according to linear 
theory S{k,z) = S{k,0)D{z)/ D{0). 

For auto-correlation, applying Limber approximation 
will change Eq. |4]to the following: 



dz 



xHz) 



l?{z)N'^{z)P{k,z) 



(6) 



For cross-correlation between two different large scale 
structure samples (be it different selection functions, red- 
shift distributions, different biases), we can write the 
cross-correlation as follows: 

Cf = I dz^--b{z)b\z)N{z)N\z)P{k,z) (7) 

where g' can have different biases, redshift dependence 
etc. 

We have not yet distinguished between the galaxy and 
the matter angular power- spectrum yet. Throughout 
this paper, we simply assume 



Cg{i) = blC(e)+N,^,,+a , 



(8) 



where Cg{t) and C(£) are the galaxy and matter angu- 
lar power spectra; bg is the linear galaxy bias, A^shot is a 



constant shot noise term which is estimated by the opti- 
mal quadratic estimator and a is a constant term that is 
fitted as a freely fioating parameter. This is a good ap- 
proximation on large scales, but breaks down on smaller 
scales; we defer a discussion of its regime of validity, as 
well as the nonlinear evoluti on o f the power spectrum to 
a later section of this paper |3.3| 

Throughout the paper, we adopt this linear redshift 
independent (within our redshift slice) bias model with 
a constant shot noise term. The bias and the shot noise 
term of galaxy sample for the various redshift slices are fit 
as extra parameters in Cosmological Monte Carlo (COS- 
MOMC; [Lewis fc Bridlej (|2002)) chains to ensure we do 
not bias our cosmological models via fixing a 
pre-computed bias. 



ig any particular 



3.2. Redshift- Space Distortions 

The position of observed galaxies can be inferred from 
their redshift, and hence the peculiar velocity along the 
line-of-sight can in principle affect our angular power- 
spectrum. So far we have neglected the effect of the pecu- 
liar velocity, i.e., the redshift-space distortion (RSD) ef- 
fect on the angular power spectrum. In the 3-dimensional 
redshift-space power spectrum measured with spectro- 
scopic surveys, the modeling of RSD is still challenging 
due to the fact that the mapping process from real to 
redshift-space is nonlinear in terms of peculiar velocit;' 
For recent efforts, se e for example Scoc cimarro' f20Q4| 
Taruya et a T ] (| 2 Q1Q[ ); [Reid fc Whit e|t2011 ); SeljaJH 
McDonald^ ( |2(J11 ). It is comparatively easy to model tHe 
RSD effect on the angular power spectrum, because the 
RSD information along the line of sight is projected out 
in the angular clustering. jPadmanabhan et al. (2007) 
formulated the RSD for the angular power spectrum at 
the linear level, and showed that the linear RSD effect 
can be seen only at large scales (^ < 20). However we 
could imagine that, if we select thin redshift slices, the 
nonlinear RSD effect may not be pro jected out and be- 
comes non-negligible at small scales. Saito et al. ( |2012 ) 
shows that such nonlinearities becomes important only 
in the case when < 0.01 at ^ > 500 but this is not the 
case here. 



We here include the linear RSD effect following Pad- 
manabhan et al. ( 2007[ ). To be co mplete, let us review 
some of the important details from Padmanabhan et al. 
(20071. 



l+^(^) = j dxN{s) [l + (5(x^,x)], 



(9) 



where we have now written the normalized selection func- 
tion as a function of redshift-space distance, s = x + v • ^ 
with the peculiar velocity component, v. Assuming the 
peculiar velocities are small compared with the thickness 
of the redshift slice, we Taylor expand the selection func- 
tion to linear order, 



dN 

iV(,)^iV(x) + ^(v-^), 



(10) 



Substituting this expression into Eq. |9) we express sep- 
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arately the 2D galaxy density field in two terms, g = 
^ g^, where g^ is the term discussed in the previous 
section, while g^ is the linear RSD correction. With the 
help of linear continuity equation, we have the Legendre 
coefficient as 



The component is given by 



(2^ 



Wl{k) . 



(11) 



where /3 is the growth parameter defined by /3 = 
d\nD/d\na/bg^ and is the derivative of the spheri- 
cal Bessel function with respect to its argument. We can 
then apply the fact that Ci = {gig})^ and calculate the 
redshift space distorted angular power-spectra. 



3.3. Non-linearities 

Non-linearities in the power-spectrum are caused by 
the non-linear evolution of components of the Universe, 
especially the late time evolution of matter and baryons. 
To capture the full extent of the non-linearities, with 
a lack of full-fledged non-linear evolution theory, one 
will need to simulate the evolution of most if not all 
of the components of the Universe. Extensive research 
and discussion have been c arried out in multiple fronts 

whether it is by perturba- 



dSanchez etl I][2QQ8 l [20091) 



tion t heory ([Carlson et aL 2009|), dark matter sin iula- 
tions (Hamaus et al.| 2010 [Heitmann et al. 2009), or 



fitting tuncti ons su^^ested by dark matter simulations 
(iSm ith et al.|2003 ) . Historically, there are a few ways to 
deal with non-linearities in utilizing power-spectrum to 
constrain cosmology, such as comparing the non-linear 
power-spectrum to the linear power- spectrum (usually 
for specific cosmological model), a nd keeping only scales 
that are believed to be linear (Tegmar k et al.|2QQ 4{|Pad- 
[manabhan et al. 2007); or utilizing the halo occupation 
model to convert a galaxy power- spectrum into a halo 
power- spectrum, which can be easily compar ed to halo 
power- spectra from dark matter simulations ( |Reid et al. 
2010); or using a variety of fitting functions developed 
Carlson et al.| (|20Q9]) to fit its observed galaxy power- 
spectra ( [Blake et al. 2010). Our project both benefits 
and suffers from the fact that it is a photometric sur- 
vey. On one hand, its BAO signal is smeared as we don't 
have accurate redshifts; on the other hand, the integra- 
tion along lines of sight ameliorates the non-linearities 
that would have been considerably stronger. Therefore, 
traditionally, angular power-spectra analysis usually only 
applies a simple cut on the an gular scale that roughly cor - 
responds to k = 0.1/cMpc~^ (Padmanabhan et al. ]|2007[ ). 
In this paper, we take a small step forward in terms of 
non-linearity treatment of the overall shape of angular 
po wer-spectrum, an d also a dopt a simila r treatment as 



3.3.1. Non-linear effects on the overall shape of the 
power- spectrum 

There is an extensive literature discussing how one can 
model the lineariti es of 3D power-spectr um over a large 
range of scales (^Sa nchez et a l. 2008; Ca rlson et al.[[2QQ9 



Hamaus et al.[[2Dl0p. 'I'his paper does not intend to ad- 



dress the issue of fully modeling the non-linearities in 3D 
power-spectrum; we do, however, take a simple model 
that happens to perform quite satisfactorily for the 2D 
angular power-spectrum. We adopt the simple linear 
redshift-independent biasing model (with shot noise sub- 
tracted for every single angular power-spectra). There- 
fore, in addition to the cosmological parameters that are 
of interest for each model, we include three extra param- 
eters for each redshift slice (6, A^shot and a) as shown in 

m 

We test the sufficiency of this model in multiple ways. 
We test this model by fitting only 2 < ^ < 150 and 2 < 
£ < 2 00 us ing simulated CM ASS mocks (as is discussed 
in Sec 3.5). We compute optimally quadratic estimated 
power-spectra of simulated data (a total of 160 realiza- 
tions from 20 independent simulation boxes, 8 lines of 
sight each), and then we compute 8 averaged (over 20 
independent simulations) power-spectra, and combined 
it with a pseudo-WMAP7 likelihood (which has the co- 
variances of WMAP7 likelihood, but with cosmologi- 
cal parameters centered on the input parameters of the 
CM ASS mocks. We find that when using above men- 
tioned model for the averaged power-spectra, in combi- 
nation with pseudo-WMAP7, we recover all input cos- 
mological parameters of the CMASS mocks for all 8 av- 
eraged power-spectra to within 1.5cr. We conclude that a 
spread over 1.5cr is reasonable. The bias parameters re- 
covered are also similar to the inpu t bias of the CMASS 
mocks as described in White et al. (2011). We therefore 
conclude that this model is accurate m recovering cos- 
mological parameters when used in the range of angular 
scales as specified above. 

We further tests this model by comparing this model 
with Hamaus et al. (2010), we found that our sim- 
ple method fits the non-linear power-spectrum derived 
from cosmological simulations quite well even up to /c = 
0.2h/Mpc. In Figure [sj we plot the non-linear power- 
spectrum from numerical simulations of halos (points 
with errorbars), while the solid lines are power-spectrum 
of various halo mass bins calculated using our simple 
model P^on-iin{k) + 4 , the model fits the non-linear 
power-spectrum quite well over a significant range in k 
even when we have not yet added the additional constant 
term a. The dashed lines show the results without the 
shot noise term for various halo mass bins. Our model of 



non-linear power-spectrum is based on HALOFIT ( [Smith 
et al. 2003), so in order to not confuse the reader, we will 



m 



Eisenstein et al.j ( [2007D ; [Blake et alT] ([2010^ for the 



non-linear treatment on the BAO scales. 



call Pnon-iin{k) by Phaiofit(^). The lower panel shows 
the ratios between a = Phh{k) — {b'^ PhaioHtik) + 1/n) and 
Phh{k) are plotted as lines. The non-linear bias is fairly 
well fit by our simple model even if we do not include the 
extra constant bias term. We decided to include the ex- 
tra constant term a to help remove the residual difference 
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Fig. 5. — To justify our choice of scale for fitting our cosmo- 
logical parameters, and the model we adopted, we show how well 
simple model such as 6^Phaiofit (^) + 1/n + a can fit fairly well 
up to /c = 0.2h/Mpc. Top panel shows the non-linear power- 
spectrum of halos in a cosmological simulations (dots with er- 
ror bars), the model being considered here 6^Phalofit(^) + 1/^ 
(solid line) and the dashed line shows what happens if we only 
use b'^{k)P^i^(k) instead. Bottom panel shows a/Phh(^) where 

= ^halofit(^) + 1/^7 we can see that the ratio is fairly consistent 
with 0.0 for a large range of mass, and it starts to deviate from 
0.0 starting k = O.lMpc/h. Therefore, we find it prudent to 
include an extra parameter a in our formalism, as we do include 
some modes at k larger than k = 0.1 Mpc/h. The different color 
lines correspond to different halo mass ranges. The largest halos 
are those with highest bias, which also gives the largest deviations 
from the model. 

between Phh{k) and 6^Phaiofit(^) + 1/^- 

3.3.2. Non-linear effects on the BAO 

We test the effect on our results of non-linear evolution 
on the smearing of the BAO feature by assuming that the 
no n-linear matter power s pectrum follows the expression 
in |Eisenstein et al.| ( [2QQ7| ): 

P(fc) =exp {-k^Y.lj2) P^^ggieik) + 

(1 - exp Pno-w^ggle{k) (13) 

where ^^i = 7.527/i-^Mpc, Pwiggie{k) is the linear 
theory power-spectrum (which includes the BAO) and 
Pno-wiggie(^) IS a smooth power spectrum, with the same 
shape as Pwiggie{k) but without any baryonic oscillations; 
w hich is computed u sin^ t he approximation described 



Eisenstein fc Hu| (1998). Both the wiggle and the 



no- wiggle part have been computed in linear theory; we 
then added to both of them the corresponding nonlinear 
ratios as a function of the scale. This approach signifi- 
cantly enhances the power in small scales. We find that 
the results are not very sensitive to the exact value of 
Sn^ prov ided that it is in the range of 5.527 to 9.527 



h~^Mpc | Eisenstein et al.|2QQ7 ). In principle, T^ni is cos 
mology dependent, and thus can change our cosmologi- 
cal constraint if it is kept as a free parameter. We have 
therefore examined our constraints on cosmological pa- 
rameters using different T^ni • We test this issue by fitting 
the full set of cosmological parameters using MCMC fit- 



FiG. 6. — We applied two different methods in calculating the 
power-spectrum (including the BAO) with non-linear effects taken 
into account, and find that it makes very little difference in the 
cosmological parameter constraints. 

ting method with COSMOMC with E^z set to 2/iMpc"^ 
higher and lower than its currently chosen value (7.527), 
and find that when we fit for a ACDM model in combi- 
nation with WMAP7, there is less than 5% change for 
any of the parameters. 

However, the addition of the nonlinear ratios is quite 
important, not only because the power in small scales 
in the angular power spectrum at high multipoles is not 
expected to be accounted for the shot noise due to finite 
number of galaxies, but also because the small shift in 
the BAO wiggles can slightly modify the best-fit shape 
of the power spectrum and hence return a different value 
of r = Qjnh. We applied two different methods in cal- 
culating the power-spectrum (including the BAO) with 
non-linear effects taken into account, and find that it 
makes essentially no difference (see Figure |6|. 

3.4. Optimal Estimation of Angular Power Spectrum 

The theory behind optimal power spectrum estimation 
is now well established, so we limit ourselves to details 
specific to this discussion, and ref er the reader to the nu- 
mero us references on the subject ( Hamilton||1997[ Seljak| 
1998" and references therein) . We also refer the reader 
to the Appendix |A| for more specific details that relate 
to our paper directly. 

We start by parametrizing the power spectrum with 
twenty step functions in (7/, 



(14) 



where the pi are the parameters that determine the power 
spectrum. We form quadratic combinations of the data. 



Qi 



1 



(15) 



where x is a vector of pixelized galaxy overdensities, C is 
the covariance matrix of the data, and is the deriva- 
tive of the covariance matrix with respect to Pi. The 
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covariance matrix requires a prior power spectrum to 
account for cosmic variance; we estimate the prior by 
computing an estimate of the power spectrum with a 
flat prior and then iterating once. We also construct the 
Fisher matrix, 



0.001 



-tr [QC-^C^-C-^] 



(16) 



The power spectrum can then be estimated, p = F~^q, 
with covariance matrix F~^. 

We also refer the reader to Appendix [A] for details more 
specific to our project. 

3.5. Tests with simulations 

To test whether the errors estimated by the quadratic 
estimator employed here are accurate or not and to test 
the results of our pipeline, we must compute the errors 
obtained via a series of simulations. 

One way to do this is to generate Gaussian random 
field using the prior power spectra for each redshift slice 
to simulate over the entire sphere. We can Poisson dis- 
tribute galaxies with probability (1 -\- S)/2 over the sur- 
vey refflon, trimmed with t he angular selection function. 
Padmanabhan et al. ( 2007 ) has tested this pipeline with 
the Gaussian random tields simulations, thus what we 
need to test here is whether the errors estimated by the 
quadratic estimator are appropriate, considering that the 
power-spectrum measurement is only minimum variance 
measurement when the field is Gaussian, which is not 
the case here. Given the non-gaussianity of the field, we 
need to determine how close we are to minimum variance 
measurement. 

As we would like to simulate our galaxy sample as 
closel y as possible , we em ployed CM ASS mock catalogs 
from|White et al. (2011) t o test the accuracy o f the op- 
timal quadratic estimator. White et al. (2011) has pro- 
duced a series of mock cata l ogs th at use the best-fit HOD 



models from White e t aL (|2011), and pop ulate a series 
of N-body simulations (jwhite et al. 201lJ). The major- 



ity of the galaxies are central galaxies living in halos of 
mass 10^^ h~^MQ. We generate 8 lines of sight from each 
corn er of each of the 20 i ndependent CM ASS simulations 
from White et al. (2011). These mock catalogs are then 
processed the sanie manner as the real data through the 
quadratic estimator code, and analyzed in the same man- 
ner as the real data set. The mock angular power-spectra 
are thus optimally estimated angular power-spectra. 

We plotted the distribution of the power- spectrum 
from each simulation that are estimated by the quadratic 
estimator code, and compare these results to the aver- 
aged error-bar of the simulation (see Figure [7|. When 
comparing the expected error to the distribution of esti- 
mated power-spectrum from each simulation and the av- 
eraged measured error from each simulation, we conclude 
that the averaged measured error is a good measure of 
the expected error. We have plotted the estimated error 
(red crosses of the middle panel of Figure [t]) by examin- 
ing the variance of the estimated power-spectrum from 
each simulation at each £-mode. We have also plotted 
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Fig. 7. — (top)The estimated (from Optimal quadratic estima- 
tor) power spectrum from 160 simulated realizations (re d crosses) 
out of 20 independent dark matter only simulations (Wh ite et al | 
2011). The green points show the averaged recovered power- 
spectrum from each of the simulations, while the blue points show 
the averaged measured error as estimated from optimal quadratic 
estimator. We conclude that the averaged measured error from the 
optimal quadratic estimator is a good measure of the expected er- 
ror, (middle) The estimated error (red crosses) is based upon the 
variance of the estimated power-spectrum from each simulation at 
each e^/-mode. We have also plotted the averaged measured error 
(green points), and it is a bit higher than the estimated error from 
the variance of the simulations (bottom panel). Nevertheless, we 
show that the quadratic estimator code can estimate the errors of 
the power-spectrum in the scale of interest here with reasonably 
high accuracy. 
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the averaged measured error (green points), and it is a 
bit higher than the estimated error from the variance of 
the simulations (bottom panel). This is probably due to 
the fact that since there are only 20 simulation boxes, 
with 8 lines-of-sights overlapping slightly in within each 
box. Therefore, the variance of simulated Ci is probably 
slightly smaller than it should be at all scales, due to the 
correlations between lines of sight. Regardless, we show 
that the quadratic estimator code can estimate the er- 
rors of the power- spectrum in the scale of interest here 
with reasonably high accuracy. It is important to note 
that at all scales of interest (to the current paper), the 
estimated error from the quadratic estimator code is not 
under-estimated. 

3.6. The Optimally Estimated Angular Power Spectrum 

The angular power-spectra estimated using the 

methodology described in Sec |3.4| are displayed in Fig- 
ure [s] In particular, we plot separately the north (Galac- 
tic), south and full angular power-spectra of these 4 red- 
shift bins (CMASS 1-4, from z = 0.45 - 0.65). We plot- 
ted the north and south separately to investigate possible 
systematic differences due to the long separation of ob- 
servation time between north and south galactic caps. 
For the scales of interest (30 < I < 150), the north and 
south are not different to prompt separate analyses. Nev- 
ertheless, this does not preclude possibility of systematic 
differences at the largest scales (at low multipole) in the 
angular power-spectrum. This is only possible, since the 
estimated power in each -^-bin is not correlated, there- 
fore a systematic difference in one -^-bin does not affect 
another. 

To test the similarity of north and south region on 
scales of interest ( 30 < ^ < 150), we find all best fit cos- 
mology parameters (with combined with WMAP7, via 
MCMC chain using COSMOMC) found by north and 
south alone respectively are consistent with each other. 
It is interesting to note that the south has smaller area 
than the north, and thus there are less information per 1- 
bin, thus the error-bars in the south is significantly larger 
than the north. It will also be discussed later in Sec [4] 
as the systematic treatment presented in this paper will 
in principle correct systematic variations even when the 
full survey is analyzed in one piece. We can also see 
the evolution of the angular power-spectra over different 
redshift slices, as it is expected. 

As shown in Figure [3) we need to investigate the poten- 
tial effects of overlapping redshift distributions. We cal- 
culate the cross-power of various redshift combinations, 
and they are shown in Figure [9j Cross-power between 
different redshift bins not only add information in terms 
of cosmology, but also from the perspective of systemat- 
ics. 

When we examine cross-power across various redshift 
bins, any difference between the measured power and 
the expected power (from galaxy auto-correlations in the 
same redshift range) can be used as a measure of the 
effects of systematics. In the top panel, there is signifi- 
cant extra power at large scale, and also negative correla- 



tions (which cannot come from galaxy auto-correlations), 
therefore, we know that there is significant systematics 
within CMASS 1. The bottom panel shows that the high 
redshift slice CMASS 4 also has substantial effects from 
systematics. 

Finally, to estimate whether it is worth including the 
cross-power of various redshift slices into the cosmologi- 
cal analysis, we performed a simple Fisher analysis. We 
calculated Fisher matrices using angular spectra from the 
four redshift bins (CMASS 1-4), with the redshift distri- 
butions given in Figure |3] A standard ACDM cosmology 
is employed to calculate the fiducial spectra. We used 
the Limber approximation (where the input power spec- 
trum was given by C AMB linear power-spectrum and 
HALOFIT) and ignored redshift space distortions. We 
employ the standard Gaussian expression for the covari- 
ance matrix of the spectra. The shot noise term was 
calculated assuming Ni = 1/n (with n being the num- 
ber of galaxies per steradian of the individual bin). Fi- 
nally, to construct the Fisher matrix, we used the range 
/ = 30—300. The parameter space is given by: 1^6, l^c, 
l^A, (78, n^, 6i, 62, 63 and {Sn refers to biases of galaxy 
sample at redshift slice N). The Fisher matrix is then 
added to WMAP7 Fisher matrix and invert to find the 
covariance matrix for the parameters. We then consider 
two cases: (1) using only the auto-spectra as observables 
and (2) using both auto- and cross-spectra as observ- 
ables. The errors on all parameters improve by less than 
5% in going from (1) to (2). We also found that ignoring 
covariances between different auto-spectra (we do include 
the covariance between auto power-spectra in the anal- 
ysis) makes less than 5% difference. This suggests that 
when we include these covariances in the MCMC, the er- 
rors will not change significantly. We therefore adopt a 
conservative approach where we don't include the cross- 
power as extra signal, but we include the bin-to-bin co- 
variance that can, in principle, be double-counted due to 
the overlap of redshift slices. 

4. POTENTIAL SAMPLE SYSTEMATICS 

Without accurately addressing known potential sys- 
tematics on the observed number density of objects in 
our sample, we cannot claim to understand its expected 
angular power-spectra, nor can we extract cosmological 
information from it. The treatment of systematics is es- 
pecially crucial for the overall shape of power-spectrum, 
since the shape does not deviate much from power-laws, 
and has no specific features such as those in BAO. The 
oscillatory nature of the BAO signal helps it from be- 
ing contaminated by any systematic signal that doesn't 
have oscillatory features. Moreover, most BAO detec- 
tion methods attempt to minimize any influence from 
the shape directly (Eisenstein et al. 2007), thus further 
shielding the BAO technique from any systematic effects. 
We will propose a novel method of systematic correc- 
tions in Sec |5] which helps mitigate the effects in power- 
spectrum caused by any systematics. 

http://camb.info/ 
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Fig. 8. — The measured angular power sp ectru m for the 4 red- 
shift bins using methodology described in Sec |3.4| We have plotted 
the full angular power-spectra, which takes into the whole sky in 
the top panel, the north (Galactic) angular power-spectra and the 
south (Galactic) angular power-spectra. Within the range of in- 
terest, the north and south angular power-spectra are consistent, 
suggesting that the systematics are at a relatively low level in the 
scales of interest, if they affect the north and south differently. 

4.1. Description of Systematics 

Here we consider not only sample systematics, but in 
particular the systematics that may contribute to extra 
(or deficit) power in the angular scale under considera- 
tion. 



Fig. 9. — The measured angular cross-correlations for the first 
3 redshift bins with the other slices. We do not show repeats of 
the cross-correlated pairs. When we examine cross-power across 
various redshift bins, any difference between the measured power 
and the expected power (from galaxy cross-correlations) can also be 
used as a measure of the effects of systematics. In the top panel, 
there is significant extra power at large scale, and also negative 
correlations (which cannot come from galaxy auto-correlations), 
therefore, there are significant systematics within CMASS 1. In 
the bottom panel, we observe the high redshift slice CMASS 4 also 
has some substantial effects from systematics at large scales. The 
CMASS 2 and CMASS 3 samples are fairly clean from systematics 
in the scales of interest. 



4.1.1. Stellar contamination and obscuration 
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Stars can in principle mimic galaxies given the right 
colors, or give rise to obscuration due to possible fore- 
ground subtractio n issues due to pre sence of a star. As it 
was pointed out in Ross et al. (2011 ) that the magnitude 
range of stars do not change its ettect on the galaxy num- 
ber density, we pick stars of magnitude 18 < r < 18.5, 
and investigate its influence on the galaxy auto power- 
spectrum. The stellar density map is plotted in Fig- 
ure Til while its auto power-spectrum is plotted in Fig- 
ure To] We use the same mask as the CM ASS samples, 
since the stars can only affect the galaxy power-spectrum 
where the two overlap. We calculate the cross-power- 
spectra between the stars and the various CM ASS red- 
shift slices, and find that there is a significant correlation 
at several angular scales (see Fig 12), mostly at large 
scales. In particular, there are strong angular correla- 
tions {£ < 10 for CM ASS 1, £ < 20 for CM ASS 2) be- 
tween stars and the galaxies at large scales, while we ob- 
serve the number of density of galaxie s is lowe r when it is 
closer to a star (as also discussed in |Ross et al.| ( |2011 )).In 
the paper, we do not include scales that are smaller than 
^ < 8 since it is much larger than the scales we are inter- 
ested in this paper. However, we will discuss more for the 
larger scales in a future publication on primordial non- 
gaussianities as the largest angular scales contain more 
information concerning primordial non-gaussianities. 

There is an extensive discussion on the stellar contam- 
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ination in the CMASS catalog in |Ross et ah] ( |201lD .The 
fundamental conclusions are that there are two separate 
effects: 1) stars can be confused as galaxies, thereby con- 
taminating the sample and inducing a positive correla- 
tion between the densities of stars and our sample and 
2) the presence of a star artificially reduces the chances 
of detecting a galaxy, th e reby i mparting a negative cor- 
relation. In Ross et al. (2011), since the band powers 
are highly correlated across bins of separating distances, 
the two effects together impart a slightly negative cor- 
relation between the number density of stars and our 
sample. In our analysis detailed in this paper, the esti- 
mated angular power-spectra are designed to have min- 
imal correlation across bins, therefore, we can see both 
positive and negative correlations over different scales {£- 
bin), as seen in Figure. 12 Given that we know stars are 
likely to contribute to trie observed number densities, we 
can take into account of the amount of contamination 
by using the above discussed techn ique. Our results are 
consistent with Ross et al. (2011) even though we do 
not detect smaller scale correlations between the stars 
and the galaxies, since the estimator employed in this 
paper produces estimates of angular band powers that 
are minimally correlated with other bins of band powers, 
while esti mates are highly co rrelated across bins in the 



analysis of Ross et al. (2011). Therefore, although the 
correlations between stars and galaxies concentrated in 
the largest scales, th ey appear in sma ller angular scale 
such as those seen inlRoss et al. (2011). 
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Fig. 10. — The angular power-spectra of various systematics we 
investigated in relation to the possible contamination to the galaxy 
power-spectra 

stars (18<r<18,5) 




Fig. 11. — The stellar density map constructed from stars of 
18 < r < 18.5. 

The sky brightness from SDSS is presented in Fig- 
ure [Tsj There are at least two scales that the sky bright- 
ness would affect the expected number of galaxies. The 
first is the width of the scan of the SDSS camera, the 
second comes from the fact that the southern cap has 
brighter sky, since we observe the south at higher zenith 
angle, thus more re-emission from the optically thin lines 
that were pumped by the sun originally. 

We use the auto power-spectrum of the sky bright- 
ness (shown in Figure 10), and its cross-correlations with 



4.1.2. Sky brightness 



galaxy densities at various redshifts (shown in Figure 14) 
to estimate the amount of contamination that can come 
from the sky. We discuss the corrections applied arising 
from sky brightness in Section [5j 

4.1.3. Seeing Variation 

Since SDSS uses a ground-based telescope at Apache 
Point Observatory, it is expected that the image qual- 
ity, primarily atmosphere seeing will affect the number 
of galaxies detected in any part of the sky. To quantify 
this, we plotted the seeing variations in the sky in Fig- 
ure [15] There is a striped pattern as different parts of 
the sky are observed in different nights, which have differ- 
ent atmosphere seeing. We use the auto power-spectrum 
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Fig. 12. — The cross-correlations between the various galaxy 
overdensity slices and stellar overdensities of 18 < r < 18.5. 

Sky Brightness (i-band) 




Fig. 14. — The cross correlations between sky brightness (i- 
band) and galaxy overdensities at various redshifts. 

seeing {i band FWHM) 




Fig. 13. — The sky brightness characterized by the i-band sky 
brightness in nano-maggies/arcsec^. 

of the seeing variations and its cross-correlations with 
the galaxy density to determine the effects of seeing on 
the galaxy overdensity clustering power. Since we can 
see that there are statistically significant but mild cross- 
correlations between the galaxies and seeing in several of 
angular band power, we correct for the seeing variations 
as discussed in Section [5] 

4.1.4. Extinction 



Fig. 15. — Seeing (image quality) map plotted in i-band Full- 
Width Half Maximum (FWHM). The stripes are probably caused 
by the changing atmosphere seeing over the observation time, as 
different parts of the sky are observed on different nights. 

We check for any residual effects on the observed 
galaxy over-densities due to Galactic extinction by com- 
puting the cross-correlations be tween the galaxy ove rden- 
sities and the extinction map ( [Schlegel et al.|p^98 ) (see 



Figure 17). Since SDSS avoids most heavily extmcted 



areas, we only have a small overlapping area where there 
is significant extinction, and galaxy data. We do not see 
a statistically significant cross-correlation between the 
galaxy (except for scale of < 20 of CMASS 1) and the 
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Fig. 16. — The cross correlations between image quality (seeing) 
and galaxy overdensities at various redshifts. 

E(B-V) (SFD map) 







Fig. 



17. — The Galactic extinction map from |Schlegel et al.| 
( |1998 |>. By comparing with the full mask of the sky, we can see that 
regions of maximum extinction, the Galactic plane is completely 
avoided. 

extinction field, therefore, we conclude that we will drop 
the galactic extinction from the list of possibly contribut- 
ing systematic effects as long as the range of interest in 
this analysis remain smaller than £ > 20. 

4.1.5. Color offsets 



Schlafly et al.| ([2Q10) reported various color offsets for 
the !SD!S!S tootp rmt, in particu l ar a n orth/south offset. 



As discussed in [Schlafly et al. (2010), the photometric 



Fig. 18. — The cross-correlations between galaxy of various red- 
shift slices with the Galactic extinction map. Since there are not 
significant correlations between the galaxies overdensities and the 
Galactic extinction map, we can drop the extinction from the list of 
potential systematic effects for CMASS 1, CMASS 3 and CMASS 
4 . The CMASS 2 sample has a significant contribution at 1 multi- 
pole, but for reasons that we will discuss in later sections, we would 
not include multipoles at ^ < 30, so this one multipole at CMASS 
2 would not affect our cosmological analysis. 

offsets can be estimated via two different ways: 1) using 
the color of stars in the imaging data; 2) using the stellar 
spectra to determines spectral classes, and then calculate 
differences between the observed and expected colors of 
stars. We adopted the later method, since it will not be 
sensitive to the intrinsic variations of stellar properties. 
However, this approach requires spectroscopy of stars, 
which is lacking in significant parts of SDSS southern 
sky. We still pursues it though, and found that there are 
no significant detection between the galaxy density map 
that overlaps with the offset map (which is lacking in 
southern sky coverage). We therefore conclude that this 
is not an important systematics in our sample. However, 
this is only a statement that at the sky which are over- 
lapped between the galaxy density map and the offset 
map, there isn't significant correlations. We will need 
more data in terms of the southern sky offsets before we 
can conclude on the effects of color offsets and galaxy 
densities. 

4.1.6. Magnitude Errors 
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Fig. 19. — The color offsets in r-i calculated from sources pro- 
vided by Schlafly et al. (2010). These are derived from the spectra 
of the stars, thus does not include an intrinsic variations due to the 
metallicity gradient of stars within our galaxy. 



As mentioned in Ross et al. (2011), the model mag- 
nitude errors in the southern cap are larger than the 
northern cap by ^ 10%, this may introduce a possible 
power excess (or deficit) at the lowest multipole. This is- 
sue however shouldn't affect any of the other multipoles 
which are the focus of our paper. 

5. NOVEL TREATMENT OF SYSTEM ATICS 

Assuming that residual systematics will exist best cat- 
alog we can construct without a serious loss of sky, how 
should we handle the remaining systematics? Without 
any evidence of possible non-linear effects of systematics 
on the observed density fields, we adopt the simplest ap- 
proach: linear relationship between the systematics and 
the observed galaxy density fields. 

We start from the following: Transforming fields from 
real space into spherical harmonic space (or 1 space in 
particular), so that < 5g5g >= Cf. 



^1 = ^1 



(17) 



where (5^ is the true galaxy density (in the £-space), and 
each 5i is the fluctuations of the map of the i-th sys- 
tematic, while ei characterizes how much i-th systematic 
contributes. With the lack of a better model, we assumes 
a linear relationship between the systematics and the ob- 
served galaxy number overdensity, but in principle this 
model can be modified to include higher order contami- 
nation due to the systematics. 

For a simple demonstration, we consider that we have 
only two systematics contributing to the observed galaxy 
density, so that i = 2. Assuming that the true galaxy 
density is unrelated to any of our systematics, we have 
the following: 

= {5l5l)+e\ {5^5,)+4 {5252)+2e^C2 {SA) (18) 

Furthermore, we have all the cross-correlations be- 
tween the systematic and the observed galaxy density 
map: 

(5°5i)=ei(5i5i)+e2(<5i<52) (19) 
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Fig. 20. — The cross-correlations between the galaxy overden- 
sities and the color offsets in g-r (top two panels) and r-i (bottom 
two panels). We can see that there is no significant correlation 
between the galaxies and the color offsets, removing these offsets 
from our list of potential systematics 
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and 



{SgS2)=€2 {6262) +ei {5i52) 



(20) 



Since we calculate all the auto- and cross- correlations 
of all the systematics (on top of the cross-correlations be- 
tween the systematics and the observed galaxy density), 
we can solve for ei and 62 (and they will be functions of 

In Figure [21] we show the cross-correlations of all 
the systematics (which contaminates the observed galaxy 
field), we find that the correlations across different sys- 
tematics are far from zero, and we must include the cross- 
correlations among systematics in our model. 

For simplicity of demonstration, we show the appli- 
cation of applying the correction from only 1 systematic 



(stars) in Figure 22 Since there is only inclusion of 1 sys 
tematic, the correction depends only on the auto-power 
of the stars and the cross-correlation between the stars 
and the observed density field. Although star is one of 
our dominant systematics, its effect on scales of interest 
{£ > 30) are quite minimal. This implies that with the 
appropriate estimator (which do not correlate powers in 
various scales), effects from systematics can be corrected 
relatively easily. This result further encourages us in 
terms of the cosmological constraining power that can 
be harnessed from future imaging surveys that will go 
deeper and wider. 

Finally, we calculate the final systematic-corrected 
power by including all three systematics which are found 
to have significant correlations with the observed galaxy 
density field. We include the auto-power (of both sys- 
tematics and galaxy density fields) and all cross-powers 
(among systematics and galaxy fields) at each angular 
scale. The final corrected power is shown in Figure [23] 

Nevertheless, as the optimal quadratic estimator pro- 
duces optimal errors and unbiased measurement only 
when the field is gaussian. In the case of highly non- 
Gaussian fields, the estimate is still unbiased, but the 
error is not optimal and may not be accurate. We can 
only test the validity of estimated error by the quadratic 
estimator if we simulate a large number of systematic 
mocks, and carry out variance tests such as the ones car- 
ried out for the galaxy density fields. 

While we understand the construction of mock galaxy 
catalogs, we do not fully understand how to construct 
mock systematic fields. We lack the 'theory' of system- 
atics fields (except probably the stellar density map). 
Therefore, it is unlikely that we can achieve optimal er- 
ror on the systematic corrections within the scope of this 
analysis. The estimated values is unbiased, but the er- 



ror can be over-estimated or maybe incorrect (Hamilton 



[1997 ). Nonetheless, when the systematic corrections are 
small, the uncertainty related to the correction cannot be 
larger than the correction itself. Therefore, we conclude 
that the most conservative way would be to include only 
power from multipoles that have relatively small correc- 
tions. Lacking a better model for the systematics, we 
adopt the following simplistic model of estimating the 
covariance of the systematic-corrected power-spectra for 
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Fig. 21. — The cross-power among systematics which are found 
to contaminate observed galaxy field. In particular, we show the 
cross-correlations between stellar contaminations (STARS), sky 
brightness (SKY) and image quality (SEEING). 



multipole bins which require small corrections. We as- 
sume Gaussianity of the fields involved, and thus using 
the following relationship: 



An 



An 



) (21) 



We modify the above equation by adding the "correc- 
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Fig. 22. — We show the corrections made by inclusion of only 
1 systematics (stellar contamination). We can see that there are 
significant corrections in large scales, but the systematic does not 
affect the smaller scales. 

tional power" due to systematics: 



fsky{2e+l) 



,-(V^C2(0 + (AQ(0)2+iVf°*) 



(22) 

This is for each = 1, so we take into account the fact 
that the Si is not 1 in ah of our bins. The quantity ACj {£) 
is the correctional power contribued by systematics. This 
method assumes the gaussianity of the fields, which is 
not a satisfactory assumption, the optimal quadratic es- 
timator can in principle project out powers that are un- 
derstood, such as those time-dependent systematics can 
be projected out in the CMB map-making. Nonetheless, 
a full modeling of the systematics and then projecting 
them out using optimal quadratic estimator is a much 
larger under-taking, which will be left to future work. 

We also show for completeness purposes the power- 
spectra of various redshift slices before and after the cor- 



rections in Figure 24 



We also compute the correlation function of our sys- 
tematic corrected power- spectr a and compare with those 
presented in Ross et clL] (|2011) and found that they are 
completely consistent with each other after systematic 
corrections. One of the systematic correla tion correc- 
tion method employed in Ross et al. (2011) follows our 



Fig. 23. — We show the systematic corrected power-spectra of 
various redshift slices when we include all three dominant system- 
atics (stars, sky and seeing). 

paper |^ and thus it is not surprising that we have 
achieved the same systematically corrected correlation 
function, even though the computation of the correlation 
function is independent. The figure [25] shows how our 
computed correlation function from our systematically- 
corrected optimally-estimated angular power-spectra is 
completely consistent (to within 1.5cr) with the measure- 
ment of Ross et al. (2011). We would also note that our 
correlation tunction shown in black lines in Figure [25] 
(our w{0) calculated using the angular power-spectra) 
has no significant large scale power, which suggests no 
exotic inflation scenario, nor residual systematics. 

6. COSMOLOGICAL PARAMETER FITTING METHOD 
AND VALIDATION 

6.1. Method 

As described in Sec |3.3[ we adopt the simple linear 
redshift-independent biasing model (with shot noise sub- 
tracted for every single angular power-spectra). There- 
fore, in addition to the cosmological parameters that are 
of interest for each model, we include three extra param- 

Even though this analysis is submitted at a later date, the 
Ross et al. (2011) paper is part of the DR8 clustering project in 
SDSS III, and thus Ross et al. (2011) has applied our method 
described here. 
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Fig. 25. — The correlation function computed from 
systematically-corrected angular power-spectra (black lines) is 
compared to the measured and systematically corrected correla- 
tion function (red squares with error bars) in Ross et al. (2011). 
The correlation function in Ross et al. 2011 also are systematically 
corrected following our method detailed in this paper. Following 
the black lines, we can see hints of the BAO at nearly all redshift 
bins, and that there is no significant large scale power seen by other 
previous analysis such as Thomas et al. (2011) prior to our DR8 
analysis. 



Fig. 24. — We show the systematic corrected power-spectra and 
the observed power-spectra (with systematics) of various redshift 
slices when we include all three dominant systematics (stars, sky 
and seeing). 

eters for each redshift shce (6, A^shot and a) as shown in 

IHl 

We have estimated that the non-Unear re dshift space 
distortion effects are minimal in our case (Saito et al. 
|2012), therefore we include the full linear redshift space 



clistortion following iPadmanabhan et al.l (|2007|) as dis 
ft 



cussed in Sec |3.2| However, calculating the full linear 
redshift space distortions requires significant time, and it 
is different from Limber approximation at / < 30; there- 
fore, we made a decision to employ Limber approxima- 
tion for multipole ranges at / > 30, and employing the 
full linear redshift space distortion calculation only at 
/ < 30. 

The measured band powers from the quadratic esti- 
mator have contributions from a range of wave numbers, 
even though they are highly concentrated in their own 
^-bin. The quadratic estimator is desi gned to compute 
ne arly a nti-correlated power spectra ( [Padmanabhan et 



al. 2003) across different multipole bins, but it still has 
a very small (<~ 5%) contribution from other multi- 
pole bins. We take this effect into account by convolv- 
ing the theory power-spectra with the window function 
before calculating the likelihood by {d — t)^C~^{d — t), 
where d represents the measured power-spectra, t repre- 
sents the theory power-spectra convolved with window 



function and C is the covariance across different bands 
and redshifts as output from the quadratic estimator. 

The fitting of the all cosmological parameters are 
do ne thr ough MCMC with COSMOMC (Lewis fc Bri-] 
die 2002). As it was discussed earlier, we do not know 
the accuracies of the error estimated using a highly non- 
Gaussian field with quadratic estimator. We know that 
it won't be biased, but the error-bar can be signifi- 
cantly mis-estimated. When the systematic corrections 
are small, we can assume that the error involved with 
the correction cannot be significant. However, at the 
largest angular scales, some of the systematic corrections 
are large enough that only proper error propagation that 
involves significant undertaking in full modeling of sys- 
tematics would provide sufficient accuracy on the error 
of the correction. There are also concerns about the va- 
lidity of the assumption of linear effects of systematics on 
the galaxy density field, especially when the corrections 
are large. We do not have indications to believe that 
the effects of systematics on galaxy density field is linear 
or non-linear. Therefore, we decided to avoid the scales 
that involve large systematic corrections, which are con- 
centrated at the large scales, thus setting the start of the 
multipole range from / = 30. 

Since our data-set is a photometric sample, the non- 
linear effects are mitigated by the fact that when we ex- 
amine the data-set, it is already integrated along all line- 
of-sights, thus decreasing its non-linearities. We there- 
fore apply the HALOFIT routine for computing the non- 
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Label 




kmax — ^max / "fiyZyy^id) 


kmax — ^max / 'fi^^mid) 






^max — 200 


^max — 150 










LRG 1 


0.475 


0.160 


0.119 


LRG 2 


0.525 


0.145 


0.110 


LRG 3 


0.575 


0.134 


0.100 


LRG 4 


0.625 


0.125 


0.094 



TABLE 2 

Descriptions of the 4 = 0.05 redshift slices, Zmid is 

THE MIDPOINT OF THE REDSHIFT INTERVAL. We ALSO SHOW THE 
CORRESPONDING kmax CORRESPONDING TO THE imax 
CONSIDERED FOR EACH REDSHIFT SLICES. 

linear power-spectrum, and limit ourselves to multipoles 
that are at relatively large scale, while non-linearities re- 
main a small effect. We choose the multipole range by 
simply testing variety of ranges using our simulations de- 
scribed in Sec 13.51 In order for MCMC chain of the sim- 
ulation to converge quickly, we average the simulations 
across different simulation boxes (so that they are not 
correlated). We tested a large range of multipole ranges: 
1 < / < 600, 20 < / < 150, 20 < / < 300, 20 < / < 200 
for example, and found 20 < / < 200 and 20 < / < 150 to 
return results within 1-sigma of the input parameters. It 
is also of interest to show the corresponding /c-limit for 
the various I ranges, since the range in k may provide 
an easier reference for the non-linearities of interest. We 
list the corresponding k range for the I range we use for 
the redshift bins considered in our analysis in Table [2j 

Combining both the low and the high / limit, we con- 
clude that 30 < / < 150 and 30 < / < 200 are both 
conservative choices for fitting of cosmological parame- 
ters. 

6.2. Cosmological Parameter Fitting Method Validation 

In this Section we present some of the tests used to 
check Crlikelihood routine for COSMOMC. To perform 
such a test we need an angular power spectrum whose 
input cosmology (i.e. the cosmological model with which 
was generated, modulo cosmic variance) is completely 
known. Therefore, the mock angular power spectrum 
described in Section 3.4, being derived from a cosmolog- 
ical simulation with initial conditions given by a known 
set of values for the cosmological parameters, provides 
an excellent testbed for our fitting routine. 

The most straightforward test is to fit each individ- 
ual angular power spectrum from the mocks and check 
that every one of them (out of the 160 available) returns 
the input cosmology. However, running 160 MCMC 
chain with only one mock power-spectra each is compu- 
tationally intensive, especially since each angular power- 
spectra has the power of ~ 0.5 — 1 actual redshift slice 
from the data, thus it will take significant time for the 
chains to converge if they converge at all. 

We therefore need to combine these mock Ci with ad- 
ditional data-set, and the mos t obvious choice be in^ the 
CMB data from WMAP7 (Larson et all [Mil). The 



WMAP7 best fit parameters however are not exactly the 
same as the simulation input parameters, thus we re- 
place the standard CMB-likelihood by a much simpler 
one, in which we compute the value of from the ac- 



tual covariance matrix from WMAP7, but not using the 
actual parameters themselves. To fully validate the fit- 
ting method, we need to use mocks that show a similar 
signal-to-noise ratio as those observed in the data. We 
combine individual mocks by using different simulations 
(and not different lines of sight in the same simulation). 

6.2.1. Building Covariance of mocks 

We compare the Gaussian covariance matrix of power 
spectrum from OQE (i.e., 'OQE covariance matrix') with 
the dispersions among 160 mocks (i.e., A/'-body covariance 
matrix'). Note that the 160 mocks are not strictly inde- 
pendent from each other, as different lines of sight share 
a small but nonzero amount of volume. To exclude the 
artificial covariance between different lines of sight, we 
derive the covariance matrix of the 20 independent mocks 
per each line-of-sight; we then average the eight covari- 
ance matrices for eight lines of sight. As a comparison, 
a straightforward dispersion among 160 mocks gives an 
almost identical result, implying that the different lines- 
of-sight share very little volume. 

In the upper panel of Figure [26j the red points are 
square roots of the diagonal elements of the OQE covari- 
ance matrix, and the black squares are from the TV-body 
covariance matrix. The diagonal elements of the OQE 
covariance matrix can be analytically calculated based 
on the smooth fit to the measured power spectrum and 
the number of independent modes for each wave num- 
ber band assuming Gaussianity, if the matrix is diago- 
nal. However, the OQE covariance matrix includes the 
effect of the window function due to the survey geometry, 
and the covariance matrix therefore is not strictly diago- 
nal and has a small anti-correlation between neighboring 
bins. Indeed, we find that there is a small deviation 
between the OQE covariance matrix and a naive Gaus- 
sian error calculation without accounting for the win- 
dow function. The difference is expected since the naive 
Gaussian error calculation does not include the effects of 
the actual survey geometry. The black dashed lines in the 
figure are the theoretical, expected error derived based 
on Gaussianity; we have rescaled it with an empirical 
boost factor of 1.1 to better match the observed disper- 
sion. The dispersions between mocks are systematically 
lower than the OQE on large scales but appears to lie 
between the OQE expectation and the boosted Gaussian 
apprixomation. 

The lower panel shows the off-diagonal elements of the 
A/'-body covariance matrix in comparison to the OQE co- 
variance, for a slice at / = 185. We observe fluctuations 
up to 20% in the measured off-diagonal terms but find no 
obvious indication that it disagrees with the OQE covari- 
ance matrix. Therefore, we conclude that the OQE co- 
variance matrix based on the Gaussian assumption does 
not underestimate the true error of the 2-dimensional 
projection of the nonlinear galaxy field. 

For the real data, we use the covariance matrix from 
the OQE for the auto-power for each redshift bin.; while 
we use the OE The upper panel of Figure [27| shows the 
square roots of the diagonal elements of the CJQE covari- 
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ance matrix of the real data (open circles) in comparison 
to the prediction based on the Gaussian prediction (after 
boosted by 1.1: solid squares) for all four redshift bins. 
The agreement is even better than the mock case, and it 
is probably due to the larger survey area of the real data 
that decreases the cross-correlation between different I 
bins that our simple Gaussian approximation cannot ac- 
cess. The lower panel shows the cross-correlation be- 
tween different I bins for four different slices of the OQE 
covariance matrix. We overplot covariance matrices for 
CMASSl (black), CMASS2 (red), CMASS3 (blue), and 
CMASS4 (magenta). Note that the covariance structure 
is identical for the four redshift bins, which is reasonable 
as the four redshift bins are subject to the same mask. 
The same structure should apply to the covariance be- 
tween different redshift bins as well. We therefore build 
the cross- covariance between different redshift bins by 
combining the diagonal elements from the Gaussian as- 
sumption and the covariance structure in the right panel 
of Figure p7l The diagonal elements are constructed us- 
ing smoothnts to the measured auto and the cross-power 
spectra of and between redshift bins and boosted by 1.1 
based on the results of the auto power spectra: 



0.001 



/sky-^mode 



Cij{l)Cij{e), (23) 



where i and j indicates a redshift slice, /sky is the fraction 
of the sky, A^mode is the number of wave modes within 
the band, and Cij{£) is a smooth fit to the auto or cross 
power spectrum; we include the shot noise contribution 
to Cij{£) in the case of the auto power spectra (i.e., i = 
j). The factor afac is the empirical factor of 1.1 that 
we introduce to match the OQE covariance matrix and 
equation [23] We use this equation to build the covariance 
between different redshift slices, while using the OQE 
covariance matrix for the covariance within the redshift 
slice. 

6.2.2. Mock test results 

With the combined average of 20 spectra each in com- 
bination with the pseudo-WMAP7, we find that the 
above model recovers all input cosmological parameters 
of the CMASS mocks for all 8 averaged power-spectra 
to within l.ba. The recovered bias parameters are also 
very similar to the in put bias o f the CMASS mocks as 
described in White et al. (2011). We therefore conclude 
that this model accurately recovers cosmological param- 
eters when used in the range of angular scale specified 
above. 

7. RESULTS 

7.1. Constraints on Cosmological Models 

The angular clustering measurement can be used to 
constrain cosmological model in several different ways: 
through standard rulers such as the matter-radiation 
turn-over scale, the baryon acoustic oscillations, or 
through large scale power which would con s train pri- 
mordial non-gaussianities (Dalai et al. 2008 Slosar et 
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Fig. 26. — Comparisons between the OQE covariance matrix 
(red circles and solid line) and the dispersions among 160 mocks 
(black squares). The black dashed lines are an theoretical, expected 
error derived based on Gaussianity after rescaled with an empirical 
boost factor of 1.1 to better match the observed dispersion. The 
lower panel shows the off-diagonal elements of the A/-body covari- 
ance matrix in comparison to the OQE covariance, for a slice at 
/ = 185. 



aL|2008). In corapanion paper (Seo 2011), we examine 
only baryon acoustic oscillations, and remove any con- 
tribution from the overall shape of the power-spectrum; 
in this paper, we include the overall shape of the power- 
spectrum and parts of the baryon acoustic oscillations 
to derive constraints on cosmological models. There is a 
companion publication on the neutrino mass constraints 
using the same angular power-spectrum (|de Putter et al. 
Ul2011|). 

Here, we choose to consider a variety of cosmological 
models, although not intending to exhaust all possibili- 
ties. We include several other data- sets to h elp break cos- 
mological degeneracies, such as WMAP7 (Larson et al. 
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Fig. 27. — The covariance matrix structure for the real data. 
The top panel compares the OQE prediction (open symbols) and 
the boosted Wick's theorem (solid symbols) for the four redshift 
bins. The two agree very well. The bottom panel shows the unifor- 
mity of the off-diagonal structure of the four redshift bins for four 
slices of the covariance matrix. The uniformity arises from the 
common mask. We therefore use Gaussian assumption and such 
uniform off-diagonal structure to build a cross-covariance between 
different redshift bins. 



2011), the "Union 2" supernova datas et (hereaft er SN), 
whicn inc ludes 557 su pernova from Hamuy et alJ(T996); 
Riess et aL (1999); Ri ess et al.|( 2QQ7); Astier et al. (2006); 

• ^(i 2006J;,Woo1i-Vasey et aL.(.2007J;.Holtzman e t. 
^clcen et al.| (2000) ; [Kessler et al.| ( |2009| , 



J ha et ar 

p05 



aT .^.^^^^ 

and Ho constraints trom usin^ 600 (Jepheids observed by 



Wide Feild Camera 3 (WFC3) published by [Riess et aL 
( 20U\ (HST). 



7.1.1. Flat CDM model with a constant equation of state 

We investigate the flat CDM model with a constant 
equation of state parameter {w) to characterize Dark En- 
ergy with the combination of our angular power-spectra 
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Fig. 28. — The 2D contour of erg and w when combining 
WMAP7 + HST + DR8. The results are consistent with w = -1. 

from SDSS-III Data Release 8 (DR8) with other data- 
sets. When we combine our DR8 observed angular 
power-spectra with WMAP7 +HST + SN dataset, we 



-1.07 ± 0.0775, 



0.2699 ± 0.0166 and 
We also combine our 



find w 

(78 = 0.85 ± 0.044 (see Fig|28|. 
systematic-corrected DR8 angular power-spectra with 
WMAP7 +HST +SN data-set, we find w = -1.064 ± 
0.0757, Qrn = 0.267 ± 0.0163. The systematically- 
corrected angular power-spectra gives consistent results 
compared with the observed angular power-spectra. 

We compare our results with other large scale struc- 
ture datasets, such as the latest large scale structure 



constraints from galaxy clustering in Blake et al. (2010), 
which has detected BAO at z ~ 0.6 usmg spectroscopic 
survey WiggleZ which includes 200,000 galaxy spectra 
over 800 deg^. They found a similar constraints on the 
equation of state of dark energy: w = —1.03 ±0.08 when 
they combined with WMAP7 + SN. This implies that 
our dataset, even though it is purely imaging data, gives 
a similar constraining power when compared to latest 
spectroscopic survey such as WiggleZ. 

We also compare our results with the BAO constraints 
from SDSS-DR7, when combined with WMAP7 + HST 

they found 



(Reid et 



3t al.||2010[ [Perci val et ~aL| [2Q1Q|) 
1.10 ± 0.14, while Montesano et a] 



w -i.iU ± U.i4, wnlTe Mont esano et al.| ( |2011[ ) used 
full shape of P(k) from SDSS-DR7, when they combined 
it with CMB + HST, they found w = -1.07±0.11. When 
we combine with the same dataset (CMB+HST), we find 
w — —1.165 ±0.12, which implies that our dataset gives a 
similar constraining power as with the full 3D DR7 spec- 
troscopic sample (at z < 0.45), while our purely imaging 
dataset is at a higher redshift range (0.45-0.65). 

7.1.2. Open ACDM model 

For an open CDM model, when combined with 
WMAP7+HST, we find = 0.00348±0.00539, improv- 
ing the accuracy over the WMAP7 +HST constraints on 
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Fig. 29. — A comparison of the constraints on the flatness 
of the Universe using our angular power-spectrum from SDSS- 
DR8 (+WMAP7 + HST) versus WMAP7+HST alone for an open 
ACDM model. 

Qk by 40% {(T{nK) = 0.0076 from WMAP7+HST); see 
Fig (29] 

When we compare our measurement to other large 
scale structure measurements such as DR7-BA0 con- 
straints; we find very similar constraints; for example, 
ohQk dis WMAP7+HST+DR7 gives Qk = -0.0023 ± 
0.0055. 

7.1.3. Flat ACDM model 

For a flat ACDM model, when combined with WMAP 
+ HST, we find Qa = 0.73±0.019, ag to be 0.817±0.023 
and Hq to be 70.5 ± 1.6 s~^Mpc~^km, which are consis- 
tent with WMAP+HST only, while improving the accu- 
racy over just WMAP+HST by ~ 5% for all parameters. 
We show the improvement on cosmological constraining 
power over the combination of WMAP7 and HST in Fig- 
ure [30l 

7.2. Companion Results 

In this paper, we are utilizing the full power spectum 
over 30 < / < 150, including both the broadband shape 
and the BAO. As a comparison, our companion paper 
Paper H, Seo et al. (2011) derives the angular diameter 
distance scale using the BAO feature over 30 < / < 300 as 
a standard ruler, while excluding nearly all the non-BAO 
information. To summarize their method and result, 
they use the angular power spectra and the covariance 
matrix shown in this paper, build a reasonable template 
power spectra based on the estimated, true galaxy distri- 
bution as a function of redshift and the concordance cos- 
mology, and fit for the distance scale, while marginalizing 
over many free parameters that account for the shape of 
the broad band. They derive DA{z)/rs = 9.2121^;^^^ at 
z = 0.54 and the result is shown to be robust against 
assumptions they make during the fitting process. Fig- 
ure |3l] summarizes the BAO fits they derived before and 
after the systematics correction. 
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In a second companion paper, de Putter et al. (2011 ), 



the angular spectra discussed in the present paper are are 



Fig. 30. — The ID marginalized constraints of Q\ and Hq when 
compared to using only WMAP7 + HST. 

used to derive a strong new upper bound on the sum of 
neutrino masses. As neutrinos suppress growth of struc- 
ture on scales above the neutrino free streaming length, 
they leave a characteris tic signature in the pow er spectra. 
To exploit this signal, de Putter et al. ( |2011 ) model the 
galaxy spectra and their neutrino mass dependence, test 
the model using mocks and show that it can be safely 
applied to the multipole range ^ = 30 — 200, while also 
considering the conservative range £ = 30 — 150. The an- 
gular clustering galaxy data are then combined with pri- 
ors from WMAP7, the HST Hubble parameter measure- 
ment, supernova distances and the (low redshift) SDSS 
BAO measurement and the resulting upper bounds are 
discussed. We quote here the conservative bound, from 
DR8+CMB+HST+BA0+SN, of Sm^ < 0.35eV at 95% 
CL, and refer the reader to the article itself for more 
details. 

8. CONCLUSION AND DISCUSSION 

We have measured the 2D clustering power spectrum of 
luminous red galaxies using the SDSS-DR8 photometric 
survey. The principal results of this analysis are summa- 
rized and discussed below. 

Using photometric redshifts, we constructed a large 
uniform sample of galaxies between redshifts z = 0.45 to 
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Fig. 31. — The BAO-only fits derived in Paper II: a stacked 
Ci/Ci^sm data of the four redshift bins before (left) and af- 
ter systematics correction(right). oiwg means the best fit 
DA{z)/rs/8.584Mpc each case. The solid red line is the best fit 
for LRG2, as a comparison, after the wavenumber is rescaled to 
z = 0.54. For more details, see Paper II. 

0.65. This probes a cosmological volume of ~ 3/i~^Gpc^, 
making this the largest cosmological volume ever used 
for a galaxy clustering measurement. The large volume 
allows us a measurement of power going from the small- 
est scales to the largest. In particular, we probe all the 
way from the smaller scales such as the Baryon Acous- 
tic Oscillations scale out to the scale of matter-radiation 
equality with one of the most accurate measurement of 
angular clustering at 2: = 0.45 — 0.65 achieved to date. 

We applied a novel approach of treating systematics 
by incorporating both the cross/auto-power of systemat- 
ics with themselves and cross-power with galaxies. This 
allows us to not only understand the impact of vari- 
ous systematics on observed galaxy number densities, 
but also allows applcation of small corrections at scales 
where the corrections are small, thus the uncertainty re- 
lated to the corrections are negligible. Since we choose 



only scales that are minimally affected by systematics, 
we expect the final cosmological constraints from both 
pre-/and post-systematic corrections are consistent with 
each other, which is indeed the case. This method can 
be improved drastically by two extra components, which 
will be left to future work. First, we should be able to 
project out the appropriate modees that are contributed 
by the systematics in similar way as it is done in CMB 
map-making (Stompor et al.|200 2|), and this can be done 
when one is estimating the optimally estimated power- 
spectra. Second, we should be able to model the dis- 
tribution of the systematics (for example, whether it is 
gaussian or not) by investigating multi-epoch data that 
is available in SDSS-DR8. In particular, this is available 
in the Stripe 82 area, which even though small (in com- 
parison to DR8 full foot-print), it contains ~ 250 square 
degrees which were multiply scanned (~ 15 — 20 repeats 
for each fields) . This will allow us to estimate the under- 
tainty of our systematic corrections properly. These two 
components will significantly improve not only our un- 
derstanding of the systematics, thus allowing us to push 
our analysis to include larger angular scales (which are 
affected by the systematics more significantly). How- 
ever, both of them requires significant undertakings in 
both data collection, code development and simulations, 
while the improvements will not dramatically change our 
cosmological interpretation in this paper, therefore, we 
leave these two components for future projects. 

For a flat ACDM model, combining our data with 
WMAP7 + HST, we find Qx = 0.7301 ± 0.019 and 
Ho to be 70.5 ± 1.6 s-^Mpc"^km. For a open ACDM 
model, when combined with WMAP7 + HST, we find 
Qk = 0.003476 ±0.00538, improved over WMAP7+HST 
alone by 40%. For a wCDM model, when combined with 
WMAP7+HST+SN, we find w = -1.071 ± 0.0775, and 
Hq to be 71.31 ± 1.65 s~^Mpc~^km, which is competi- 
tive with the latest large scale structure constraints from 
spectroscop ic surveys such as those by WiggleZ ( [Blake] 
et al.|20l"Q |, and SDSS DR7 spectroscopic surveys, espe- 



cially in analys is led by ([Reid et al.||2010[ Percival et al. 

al. 201ip . This results implies that 



|201 0| FMontes ano et 



our dataset, even though it is purely imaging data, pos- 
sesses a similar constraining power as the spectroscopic 
surveys such as WiggleZ or SDSS-DR7. What we lack in 
redshift precision, we compensate by shear volume. This 
suggests that future and upcoming imaging surveys such 
as PanStarrs0 DES^and LSSTRcan achieve signif- 
icant cosmological constraints via large scale structure 
clustering even when compared to other spectroscopic 
surveys. 

This is the Paper I of the project, which mostly de- 
scribes the construction of the data- set, treatment of sys- 
tematics, estimation of the angular power-spectra and fi- 
nally using the overall shape of the angular power-spectra 
over a large range of angular scale to derive constraints 
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on our cosmological models. We refer readers to Paper 
II (|Seoll2Qlll) of the project, which uses only the Baryon 
Acoustic Oscillation feature to fit for various cosmolog- 
ical parame ters. We also refer to Paper III (de Putter 
et al. |[2Q11 ) of the project, which uses the overall shape 
ot the power-spectra to fit for various neutrino models. 
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APPENDIX 
QUADRATIC ESTIMATOR 
Consider a Gaussian random field Xi with (xi) =0 and covariance 



N 



(a) 



^=1 



We wish to form an estimator, p^, of which is quadratic in the data 

Pa ^ ^ ^ij '^i'^j 



(Al) 



(A2) 



where Q is symmetric. Requiring the estimator to be unbiased 



(Pa> = Pc. 

For Gaussian Xj the covariance of pa is 



tr 



6"" , 6a = tr 



Q(«)C(0) 



Gov = J2Q\fQi^i^ [CikCji + CuCjk] = 2tr [cQ^^'^CQ^^^ 

ijkl 



(A3) 



(A4) 



This problem is easiest if we consider a single parameter at a time, with all other parameters held fixed (and absorbed 
into C^^^ ). Thus we wish to minimize 



Taking derivatives with respect to the components of Q^"^ gives 
or 



(A5) 

(A6) 
(A7) 
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If the dependence of C on is not linear then we can use a Newton-Raphson iteration where now C^"^ is the derivative 
of C evaluated at the current best value of p. Iterating, by replacing p ^ p -\- Sp until the best-fit Sp = 0, results in a 
maximum likelihood solution. In practice, it only takes a few iterations to achieve the maximum likelihood solution. 
This approach also results in another fact that is under-appreciated in t he lit erature. The above choice of Q (which 



can have a slightly different form, see Table 4 of iPadmanabhan et al. (2003) for more details) produce error-bars 



that are anti-correlated across different band powers. In this paper, we include the window function (which is mostly 
affected by the mask) before we compare the observed power and the theoretical power. 



